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Abstract 

We propose explicit forms for the steady state distributions governing fully developed turbulence 
in two and three spatial dimensions. We base our proposals on the crucial importance of the area 
and volume preserving diffeomorphisms in the space of velocities. We argue that these distributions 
can lead to the relevant (Kolmogorov and Kraichnan) scaling laws. 
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In this note we propose universal steady state distributions for fully developed turbulent 
flows in two and three dimensions (2d and 3d) building on the recent work which aims 
to relate modern quantum field theory, string theory and fluid dynamics [l). The basic 
dynamical equation for turbulent flow is the non-linear Navier-Stokes equation (we use the 
sum convention throughout the paper) 

p{d t Vi + VjdjVi) = -dip + vd?Vi, (1) 

with the incompressibility condition diVi = 0, where vu is a component of the velocity field 
of the flow, p is pressure, and p is the fluid density |2j. In this note we are interested in 
fully developed turbulence, or turbulence in the limit of infinite Reynolds number R. As 
R = Lv/v goes to infinity, with v = y/ViVi, whereas L is a characteristic scale and v is 
the kinematic viscosity, we effectively have the limit of vanishing viscosity v — > 0. In this 
regime, all the various possible symmetries are restored in a statistical sense, calling for a 
probabilistic description of what is in essence a deterministic system (strongly dependent 
on the boundary conditions). Therefore, in computing correlators of the velocity field, we 
should use the statistical and quantum field theoretic descriptions of turbulence [3]. 

In particular we are interested in the form of a steady state turbulent distribution P(v) 
which should imply the various famous scaling laws (the Komogorov scaling in 3d |4J and 
both Kolmogorov and Kraichnan scaling laws in 2d [5]). The two-point correlators deter- 
mined by these scalings should be computable from P(v) as follows 

_ j Dv Vi {l)v^)P{v) 
M0«i(0)> = — — • (2) 

The fundamental question is: what is P(v)7 In this note we wish to argue that the answer 
to this question is: 

P(v) = exp[-S K (v)}, (3) 

where we assume that the expression for the Sk{v) is local and universal (see below). Both 
of these assumptions might be challenged: we can a priori expect non-local factors in Sk{v). 
Also, the relative locations of the Kolmogorov and Kraichnan distributions in 2d turbulence 
energy spectra depend on the (location of the) forcing, which challenges universality. We 
also note that the notion of universality might be challenged by an a priori dependence 
on the boundary conditions. We will simply assume that universal distributions exist in 
the rest of this note, in spite of these caveats. Finally, we note that if universality applies 



to turbulence, it is expected at short distance [3J, as opposed to the usual long distance 
universality associated with quantum field theory and critical phenomena from equilibrium 
physics. 

The central message of this note is that Sk is what we call the Kolmogorov distribution 
in 3d, determined by the effective action for a 3d gauge theory based on volume preserv- 
ing diffeomorphisms. Alternatively, Sk is what we call the Kraichnan distribution in 2d, 
determined by the effective action for a 2d gauge theory based on area preserving diffeomor- 
phisms. In both cases the diffeomorphisms act in velocity space. In both cases we crucially 
use the Galilean symmetry in the space of velocities, so that the effective respective actions 
involve only derivatives in v^. 

We pause briefly to remark that in order to fully characterize systems far from equilibrium 
knowledge of the probability currents is needed in addition to the steady state probability 
distributions 6]. However, a discussion of these currents is outside of the scope of the present 
paper. 

Even though the puzzle of fully developed turbulence is essentially a strongly coupled 
problem, we motivate our discussion (and ultimately, our proposals) by some rather ele- 
mentary observations at weak coupling. We start with 2d and then move on to 3d. In 2d 
the Lagrangian description of the fluid is generated by the following Lagrangian [7| (in the 
notation of [3]) 

L2 = jd'yp [^-V( Po \^-\)]. (4) 

where po is the constant density in the co-moving coordinates and the real space density is 
Po i fx I ' wnere I ff I denotes the Jacobian connecting the co- moving coordinates y and the con- 
tinuum fields Xi(y, t). The Lagrangian L 2 is invariant under area preserving diffeomorphisms 
in 2d fl: 

y'i = yi + fi(y) ; Sx a = —^fi{y), (5) 

oyi 

where for area preserving diffeomorphisms (e^- being the Levi-Civita symbol in 2d) 

and where A(y) is an arbitrary function, which generates these "gauge" transformations. 
This equation in turn leads to the 2d Poisson bracket action: 

dx a dA(y) _ ,. . ( 
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Suppose we consider small motions of the 2d fluid 

Xi = yi + eeijAj, (8) 

where e is the small coupling (and where we have absorbed the factor of po compared to 
To linear order the area preserving diffeomorphisms become the usual gauge transformations 

SAi = <9,A, (9) 

and the quadratic Lagrangian that is invariant under this transformation is just the usual 
Maxwell Lagrangian 

d2 vlA 2 - (v x If). (io) 

The inclusion of non-linear terms can be done by extending the linear gauge transforma- 
tions to their non-Abelian completion. This leads to a non-Abelian gauge theory, with the 
following map between the full area preserving diffeomorphism group and the non-Abelian 
transformations generated by a commutator of two matrices 

{A,x}^[\,X\. (11) 

Here we have used the standard [8] mapping between the Poisson brackets and commutators 
of infinite square matrices A and X. The corresponding gauge theory action reads as 

sa= ^I d2y[A * 2 ~ ( 9iAj ~ 9jAi + [Au Aj])2] - (i2) 

We will return to this gauge theory action in what follows. In 3d we can extend the presen- 
tation of [7] by considering 

Ls = jd 3 yp [^-V( Po \^\)}. (13) 

The Lagrangian L 3 is invariant under volume preserving diffeomorphisms 

dx a dA 1 (y)dA 2 (y) f . 

ox a = -^—e ijk — . (14) 

dy { dy,j dy k 

where Ai and A 2 are the generators of the volume preserving gauge transformations. This 
is equivalent to the following Nambu bracket [9]: 

8x a = {Ai, A 2 ,x a }, (15) 



4 



where, by definition 

{A, B, C} = e abc d a Ad b Bd c C. (16) 
Here A, B, C are three functions of three spatial coordinates x,y,z. This classical bracket 



seems to be naturally generalized to a triple algebraic structure [9- 

[Ai , Aj , A%\ = e a h c A a AhA c (17) 

Suppose we consider small motions of the 3d fluid in analogy with the 2d case: 

Xi = Hi + feijkBjk, (18) 

where B^ = —B^j and / denotes a small coupling. The "quantization" of the volume 
preserving diffeomorphisms is a more involved problem j^, [h|. Nevertheless, even in this 
case we expect a mapping between the full volume preserving diffeomorphism group and the 
3-bracket 

{A 1; A 2 ,a;}^ [A 1; A 2 ,X], (19) 



where Ai and A 2 are the appropriate matrix realizations of Ai and A 2 10]. The 3d action 
invariant under the linear part of this transformation is 

Sb = ^pJ d3y ^ j ~ ^ dlBjk + djBkl + dkB ^- (20) 

Note that the explicit linear volume preserving transformations are 

5B rj = d i A 1 d j A 2 - diAzdjAu (21) 

where Bij is dual to a 3-vector in three dimensions, = \eijkBjk- 

The above linear analyses seem removed from such a strongly coupled problem as fully 
developed turbulence [l^ . Still the linear analysis is useful, because in the stationary case we 
can formally replace Xi — > V{ and talk about velocity Lagrangians. That the corresponding 
Lagrangians should be given in terms of the derivatives of velocity is fixed by the Galilean 
symmetry, t>j — > i>j + Ui, for a fixed velocity with component Uj. Thus, if we follow the linear 
analysis in the velocity space and drop the time dependence because of our interest in the 
stationary distributions, in 2d we have the following natural theory that is consistent with 
area preserving diffeomorphisms involving vf. 

4 2) = ^ / d 2 v[( 9 ^ ~ d W + K (22) 



where g 2 denotes the appropriate dimensionful coupling constant and i,j = 1,2 (the factor 
is canonical). This theory is invariant under 

6vi = {F,Vi}, (23) 

where F generates area preserving gauge transformations in V{ space. We want to argue 
that this is the stationary distribution we have been looking for in 2d. First, we attempt 

to justify this guess on more general grounds: Obviously we have the Galilean invariance 

(2) 

in the inertial range v j — > Vi + u% (where Ui is constant) because S K is a functional of the 

(2) 

derivatives of v*. Second, it is reasonable to expect that S K is governed by the conserved 
quantities. (Recall the case of equilibrium statistical mechanics, i.e. the Boltzmann-Gibbs 
distribution, where S is simply the energy, an additive conserved quantity). In our situation 
we have vorticitv (Co = V x v) squared and velocity squared, as the natural conserved 

quantities (2), 0, a]. Therefore, if we start with the vorticity squared (u 2 ) term we see that 

(2) 

this is really the quadratic part of our guess for S K . This term is invariant under the linear 
gauge transformations in the space of velocities. However, by going to real space, we may 
invoke the full non-linear group of area preserving coordinate transformations generated by 
the Jacobian (the Poisson bracket) in 2d space. Then by concentrating on the steady state 
regime we may claim the same symmetry in the v space which would lead us to the above 
proposal. 

This action should be compared to the 2d Yang-Mills theory action Sa given in equation 
( fl2|) . We can use the covariant derivative di — > Di = di + A i} to rewrite the usual 2d Yang- 
Mills theory Lagrangian as [Di,Dj] 2 = (diAj — djAi + [Aj, A,-]) 2 . This procedure can be 
immediately generalized to 3d so that we have a theory consistent with volume preserving 
diffeomorphisms in the space of Vi. By using the covariant derivative di — > D\ = di + v «, we 
can immediately rewrite the above guess for the 2d action 

S k=^1 J dMW.-DJ} 2 ], (24) 

and then extrapolate our 2d proposal to the natural proposal for the Kolmogorov distribution 
in 3d 

s « ) = w I d3y[{D ^ D ^ DV ^- (25) 

Here g 2 denotes the appropriate dimensionful coupling constant. The crucial non-linear part 
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which replaces J d 2 y({vi,Vj} 2 ) is 

This action is fixed now by the volume preserving diffeomorphisms in the 3d velocity space 

^ = {^1,^2,^}, (27) 

where F\ and Fi generate the volume preserving gauge transformations. The argument 
which leads to this proposal is just the repetition of the argument we have presented for the 
2d Kraichnan distribution. 

These are thus our explicit proposals for the turbulent distributions in 2d and 3d: they 
are encoded in the expressions for and S^. Do these educated guesses give the correct 
results? We concentrate on the case of 2d turbulence. Let us remember that for the Kraich- 
nan scaling in 2d we want to derive (vi(l)vj(O)) ~ l 2 S{j, and for the Kolmogorov scaling law 

2 

(vi(l)vj(0)) ~ l^Sij. In order to discuss the validity of our proposal, we turn to the natural 
loop variables introduced by Migdal [13J]. The natural loop variable, the Migdal loop [kJ 

W m (C) = (exp(-- f VidXi)), (28) 
v Jc 

(where C is a contour and the viscosity v plays the role of an effective H Q) allows us 
to rewrite the Naiver-Stokes equations [13j as an effective Schrodinger equation 

iudtW(C) = H C W{C), (29) 



with the appropriate loop equation Hamiltonian He [131 ] . In 3d Migdal observed a self- 
consistent scaling solution of this equation in the v — > limit (a WKB limit in this problem) 
Wm(C) ~ exp(— [-^]f), which precisely corresponds to the Kolmogorov scaling. In 2d, the 
Kraichnan scaling leads to the area law [ij for the Migdal loop Wm(C) ~ exp(— A/A ). In 
the above gauge theory of 2d velocities this area law is very natural, because of the fact 
that in the corresponding 2d Yang-Mills theory, the Wilson loop (the natural analog of the 
Migdal loop) 

W{C) ee (exp(- / Aidxi)), (30) 

n 

obeys the same area law p4 W{C) ~ exp(— A/Aq). Given the precise structural mapping 

(2) 

between the 2d Kraichnan theory (defined by S K ) and the 2d gauge theory, this area law 
scaling for the loop variables should be obeyed, and the Kraichnan scaling should thus follow 



from our proposed steady state distribution. This is one of the central observations of this 
note. 

To summarize: Given the above analogy of the proposal for the Kraichnan distribution 
and what we know about the 2d Yang-Mills theory, Kraichnan's scaling implies the area law 
for the Migdal loop operator. In view of this, and assuming locality and universality, our 
proposal for the Kraichnan distribution, viewed as the gauge theory of 2d area preserving 
diffeomorphisms in velocity space, is essentially unique. 

How about the Kolmogorov scaling in 2d? Does this scaling also follows from our proposal 
for the universal steady state distribution of 2d turbulence? Here we can only offer a very 
heuristic argument. Preliminary numerical investigations [ijj] seem to indicate that the 
Kolmogorov scaling does agree with our proposal. 

We start our intuitive argument regarding the Kolmogorov scaling in 2d by noting that 
in the small momentum regime the linear term in our 2d gauge theory in velocity space 
dominates and this is where we expect Kraichnan's scaling. In the high momentum regime 
the Jacobian term (the Poisson bracket) dominates. In some sense this is an effective "mass 
term" (or the restoring force quadratic term from the harmonic motion) that we could 
associate with the kinetic energy and the velocity squared term. Given the fact that the 
energy cascade leads to the Kolmogorov scaling |2j we should thus expect the Kolmogorov 
scaling at high momenta. Thus, naively we expect Kraichnan's scaling at low momenta 
and Kolmogorov's scaling at high momenta 16j. (The issue of the possible dependence on 
forcing in the mid momentum range makes this intuition in principle more complicated j^j.) 
Nevertheless, by continuing with this picture, we could envision a large Migdal loop and 
apply the area law there (in the large distance, low momentum limit). In the center of this 
large loop we envision a much smaller loop, separated by some radial scale from the large 
one. The area in between the large and the small loops sweeps an effective volume as the 
large loop sweeps its area. We can then dimensionally translate the area A of the surface 
swept by the large loop into this effective volume V, A = V 2 ! 3 . This effective volume is in 
turn the area of the smaller loop times the distance between the scales. When we shrink 
this radial distance to one, we get the area of the smaller loop (which we might call the 
Kolmogorov loop) to the power 2/3, i.e. A 2 ^ 3 , which would correspond to Kolmogorov's 
scaling of the Migdal loop (i.e. exp(— (A/A ) 2 ^ 3 )). Needless to say, this is a very heuristic 
picture. The big caveat in this obviously qualitative argument is the issue of forcing in the 
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mid momentum range. This intuitive argument has been tested in the preliminary numerical 



simulations 



151 ] . The results are qualitatively encouraging and will be discussed elsewhere 



There exists a natural reduction from 3d to 2d when considering the reduction from 
volume preserving diffeomorphisms to area preserving ones 

{Ai, A 2 , x a } ->■ {Ai, x a }, (31) 

(provided A 2 = 1) and the corresponding reduction of the 3-bracket to the 2-bracket 

[\ u \ 2 ,X]^[\ h X], (32) 

(provided A 2 = 1). Given the above qualitative argument concerning the emergence of the 
Kolmogorov scaling (via the exp(— (A/A ) 2 ^ 3 ) scaling of the Migdal loop) and this reduction 
procedure we would expect that our proposal for the 3d steady state distribution leads imme- 
diately to the Kolmogorov scaling. Of course, only numerical tests (which are considerably 
more difficult in 3d as compared to 2d) could confirm this expectation. Given this reduction 
we can also propose an interpolation between the 2d and 3d steady state distributions for 
quasi-2d fluids of finite thickness by looking at 

4 2 ' 3) = ^- 2 j d 2 y[{d iVj - d jVi + {v hVj }f\ + ±J d?y[{v hVj ,v k }) 2 }. (33) 

Here one should perform the expansion of the 3-bracket in the third anisotropic direction 



{A, B, C} = A{B, C} + B{C, A} + C{A, B}, (34) 

where the dot derivative is with respect to the anisotropic third direction. 

Finally, we note another general caveat: Presumably turbulent laws are generated by 
some boundary conditions (or by stirring) that are responsible for the cascades and the 
"energy anomaly" in the Kolmogorov case 0, ^]. How do we see this in our proposal? 
We should remember that we need to compute the functional integral over the velocities 
(J Dvexp(-SK)) and the measure in that functional integral might transform anomalously 
once we start worrying about the boundary conditions and not only the symmetries. That 
is, the velocity space area preserving diffeomorphisms in 2d and volume preserving diffeo- 
morphisms in 3d. We expect that our discussion is correct as long as boundary effects are 
neglected. 



In conclusion, we have proposed universal statistical distributions for fully developed 
turbulence in 2d and 3d. The proposed steady state distributions for turbulence in 2d and 
3d are essentially fixed by the Galilean invariance and the area and volume preserving gauge 
transformations, in 2d and 3d, respectively. We have argued that the Kraichnan scaling 
in 2d should follow based on what is known about the non-Abelian gauge theory in 2d. 
We have also argued that the Kolmogorov scaling should be consistent with our proposal, 
based on a very heuristic picture. Obviously, only a thorough numerical study can show 
whether our proposals are true or not. We will address these numerical issues in an upcoming 
publication. 
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